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Abstract 

We report on the experimental study of an exceptional point (EP) in a dissipative microwave 
billiard with induced time-reversal invariance (T) violation. The associated two-state Hamiltonian 
is non-Hermitian and non-symmetric. It is determined experimentally on a narrow grid in a 
parameter plane around the EP. At the EP the size of T violation is given by the relative phase 
of the eigenvector components. The eigenvectors are adiabatically transported around the EP, 
whereupon they gather geometric phases and in addition geometric amplitudes different from unity. 

PACS numbers: 02.10.Yn,03.65.Vf,11.30.Er 
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We present experimental studies of two nearly degenerate eigenmodes in a dissipative 
microwave cavity with induced Tviolation. Due to the dissipative nature of the system 



the associated Hamiltonian is not Hermitian 



[i 



4] and thus may possess an exceptional 



point (EP), where two or more eigenvalues and also the associated eigenvectors coalesce. 
In contrast, at a degeneracy of a Hermiti 
the eigenvectors are linearly independent 



In contrast, at a degeneracy of a Hermitian Hamiltonian, a so-called diabolical point (DP), 

5, Is]. The occurrence of exceptional points 2], 7] 
in the spectrum of a dissipative system was studied in quantum physics Q] as well as in 
classical physics It has been demonstrated, that this is not only a mathematical but also 
a physical phenomenon. The first experimental evidence of EPs came from flat microwave 



cavities 
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13] . which are analogues of quantum billiards [14] . Subsequently EPs were 



observed in coupled electronic circuits [15[ and recently in chaotic microcavities and atom- 
cavity quantum composites jl^]. The present contribution is the first experimental study of 
an EP under Tviolation. It is induced via magnetization of a ferrite inside the cavity by an 
external field B 17| . Tviolation caused by B is commonly distinguished from dissipation 
181 ] . For a non-dissipative system with broken Tinvariance (B 7^ 0) the Hamiltonian is 
Hermitian. For a dissipative system the energy is not conserved such that for B = it is 
described by a complex symmetric Hamiltonian. This case is referred to as the T-invariant 
one [18 ]. 



To realize a coalescence of a doublet of eigenmodes in the experiment two parameters are 
varied. Since the doublets considered are well separated from neighboring resonances, the 
effective Hamiltonian is two-dimensional. Its four complex elements are determined on a 
narrow grid in the parameter plane, thus yielding an unprecedented set of data. This allows 
to quantify the size of Tviolation (ititomeasure to a high precision the geometric phase 
6[ and the geometric amplitude [Jj, |19|, |20] that the eigenvectors gather when encircling 
an EP. An earlier rule H| on the encircling is generalized to the case of Tviolation [l|. 

The experimental setup is similar to that used in 
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121 ] . The resonator is constructed 



from three 5 mm thick copper plates, which are sandwiched on top of each other. The center 
plate has a hole with the shape of two half circles of 250 mm in diameter, which are separated 
by a 10 mm bar of copper except for an opening of 80 mm length (see the inset of Fig. [I]). 
The opening allows a coupling of the electric field modes excited in each half circle, which 
is varied with a movable gate of copper of length 80 mm and width 3 mm inserted through 
a small opening in the top plate and operated by a micrometer stepper motor. The bottom 



of this gate is tilted to allow a precise setting of small couplings. When closing the gate 
it eventually plunges into a notch in the bottom plate. Its lifting s defines one parameter 
with the range mm (no coupling) < s < 9 mm (full coupling). The left half of the circular 
cavity in Fig. [TJ contains a movable semicircular Teflon piece with diameter 60 mm and 
height 5 mm connected to the outside by a thin snell operated by another stepper motor. 
Its displacement with respect to the center of the cavity defines the other parameter 5. 

Two pointlike dipole antennas intrude into each part of the cavity. A vectorial network 
analyzer (VNA) Agilent PNA 5230A couples microwave power into the cavity through one 
antenna a and measures amplitude and phase of the signal received at the same (reflection 
measurement) or the other (transmission measurement) antenna b relative to the input 
signal. In this way the complex element Sb a of the scattering matrix S is determined. To 
induce T violation, a ferrite with the shape of a cylinder of diameter 4 mm and height 
5 mm is placed in the right part of the cavity and magnetized from the outside by two 
permanent magnets. They are mounted to a screw thread mechanism above and below the 
resonator and magnetic field strengths mT < B < 90 mT are obtained by varying their 
distance. To automatically scan the parameter space spanned by (s,5), the two stepper 
motors and the VNA are controlled by a PC. The four S"-matrix elements Sf, a (f), {a, b} e 
{1,2} are measured in the parameter plane with the resolution As = A5 = 0.01 mm 
and a frequency step Af = 10 kHz. The frequency range of 40 MHz is determined by 
the spread of the resonance doublet. Lack of reciprocity, i.e. S12 7^ S21, is the signature 
for Tviolation l?]]. Figure [1] shows two typical reflection spectra, i.e. \S aa \ 2 with a = 1,2. 
Additional measurements include neighboring resonances, which are situated about 250 MHz 
away from the region of interest, in order to account for their residual influence. In the 
considered frequency range the electric field vector is perpendicular to the top and bottom 
plates. Then, the Helmholtz equation is identical to the Schrodinger equation for a quantum 



billiard of corresponding shape Jjj, |2l| . Consequently, the results provide insight into the 
associated quantum problem. 



In previous experiments III, Il2j an EP was located by determining for each parameter 
setting the real and imaginary parts of the eigenvalues as the frequencies /i^ and the widths 
of the resonances from a fit of a Breit-Wigner function. This procedure, however, fails 
at the EP because there the line shape is not a first order pole of the S matrix but rather the 
sum of a first and a second order pole 0, Q. Therefore, the coalescence of the eigenvectors 



should be incorporated in the search for the EP [ljl, For this we determine the two- 
state Hamiltonian H and its eigenvalues and eigenvectors explicitely for every setting of the 
parameters s,8, B from the measured ^-matrix elements via the method presented in 
There, we showed that for a resonator with two pointlike antennas a resonance doublet is 
well described by the two-channel S matrix S(f) = 1 — 2-niW\] 1 — H)~ l W . The matrix 
W = (Wfxa) couples the resonant states /x = 1,2 to the antennas a = 1,2. It is real since 
that coupling conserves T . The Hamiltonian H comprises dissipation in the walls of the 
resonator and the ferrite and T violation and thus is neither Hermitian nor symmetric, which 
implies 7^ S21. Its general form is given as 



H=\ 61 *&-«^|. (1) 

#f 2 + iH& e 2 



The quantities e± ± e 2 , Hf 2 and H^ 2 are complex expansion coefficients with respect to the 
unit and the Pauli matrices. The ansatz for S was tested thoroughly by fitting it to S- 



matrix elements measured with and without magnetization of the ferrite. In the latter case 
the antisymmetric part H^ 2 vanishes and H coincides with that used in HQ. Fitting 
the S matrix to the four measured excitation functions S& a (/) yields the matrices H and 
W up to common real orthogonal basis transformations. We define the basis such that 
{Hf 2 + iH^ 2 ) I {Hf 2 — iH^ 2 ) = exp(2ir) is a phase factor. Thus, T violation is expressed by a 



real phase r. This is usual practice in physics, e.g., for nuclear reactions [22[, and for weak 



and electromagnetic decay [23 1. 



The eigenvalues of H in Eq. ([T]) coalesce to an EP when Hf 2 + H^ 2 + [e\ — e2) 2 /4 = 
but not all three terms equal to zero. Figure [2] shows one example for a search of the 
EP. Keeping s, B fixed at s — 1.66 mm and B = 53 mT, the real and imaginary parts of 
the eigenvalues, fj and Tj are shown as functions of 5. At 5 = 41.25 mm the encounter 
of the eigenvalues is closest. To determine the parameter values for the coalescence of 
the eigenvectors jr^), k = 1,2, their components r^r^ 2 and the ratios = r^i/r^ = 
|^fc|e** fc are determined [6j. The right part of Fig. |5] demonstrates that the encounter of the 
eigenvectors is also closest at 5 = 41.25 mm. Indeed, as will be evidenced in further tests 
below, both the eigenvalues and the eigenvectors cross at these parameter values. The EP 
is located at (sep,8ep) = (1-66 ± 0.01,41.25 ± 0.01) mm. At the EP the only eigenvector 



of if is 




The ratio of its components is a phase factor. For T-conserving systems the phase is $ep = 
7r/2 [12] as confirmed by Fig. [3] at B = 0. With T violation the phase equals $ep = t + n/2 
and thus provides a measure for its size r. Fi gur e [3] shows that with increasing B the 



parameter r goes through a maximum. As in Ref. (17| we identify this with the ferromagnetic 
resonance resulting from the coupling of the rf magnetic field to the spins in the ferrite. The 
T-violating matrix element iH^ 2 has been expressed by a resonance formula which yields 
the solid curve for $ep(B). 

Panel (a) of Fig. H] shows in a neighborhood of the EP at B = 53 mT the differences of the 
complex eigenvalues, f\$ (blue, left to EP) and (orange, right to EP), panel (b) those of 
the phases $12 (orange, left to EP) and of the moduli (green, right to EP). The darker 
the colour the smaller is the respective difference. The darkest colour visualizes the curve 
(branch cut) along which it vanishes. The left and right panels of Fig. [2] are cuts through 
the respective panels of Fig. 2] at s = 1.66 mm. We observe that |/i — /2I and |$i — $2! are 
small and thus visible only to the left of the EP, whereas \Ti — T 2 1 and \ \i>i\ — | z^ 2 1 1 are visible 



only to its right. Thus the branch cuts all extend 



directions. This proves that this point is an EP 



rom one common point into opposite 



16|- 



For systems with T violation the geometric phase 7 gathered by the eigenvectors around 
an EP is predicted to be complex yielding a geometric amplitude e _Im7 7^ 1 [l]. To check 
this we choose contours around the EP for the six values of the magnetic field B consid- 



ered in Fig. [3J One example is shown in Fig. HI As proposed in 20| it consists of two 



different 



oops. The path is parametrized by a real variable t with initial value t = 0. In 



24j and 11( the geometric phase gathered around a DP, respectively an EP, was obtained 
for just a few parameter settings, because the procedure - the measurement of the elec- 
tric field intensity distribution - is very time consuming. We now have the possibility to 
determine the left and right eigenvectors, (lj(t)\ and \rj(t)), j = 1,2 of ff in Eq. (0Q) on 
a much narrower grid of the parameter plane. At each point t of the contour they are 
biorthonormalized, such that (lj{t)\rj(t)) = 1. Defining in analogy to the T-conserving case 



11 1 B = ^^I^Hf* + Hf 2 2 , tan 6 = a/1 + B 2 - B yields for the right eigenvectors fl 

, e- ir / 2 cos#\ ( -e- iT / 2 sm9 , 

\n(t)) =\ . , \r 2 (t)) = . (3) 

e JT / 2 sin# / \ e iT / 2 cos# 



As in JjJ the EPs are located at 1 + B 2 = and encircling the EP once changes 9 to #±7r/2. 
The T-violating parameter r is not constant along the contour, even though the magnetic 
field is fixed. In fact, it varies with the opening s between both parts of the resonator, 
because the ferrite is positioned in one of them, see Fig. [TJ The position of the EP does 
not depend on the value of r. Thus, the space curve (s(t), 5(t), r(t)) winds around the line 



(sep, $ep, T {t)), see Ref. |l9|. Since r has no singular points in the considered parameter 
plane it returns to its initial value after each encircling of the EP. As a consequence, the 
eigenvectors \ri^{t)) follow the same transformation scheme as in the T-conserving case. 
After completing the first loop at t — t\ we have |r 1 (t 1 )) = |r 2 (0)), |r 2 (£i)) = — |rj(0)) and 
after the second one at t = t 2 we find |t 1;2 (^2) ) = — 1^1,2(0)) . 

The biorthonormality defines the eigenvectors and \rj) up to a geometric factor, so 
that (Lj(t)\ = (lj(t)\e-^ {t) and \Rj(t)) = \rjt)) e^ (i) . The geometric phases jj(t) are 
fixed by the condition of parallel transport |5|, |25j], (Lj(t)\-^Rj(t)} = 0. This yields 
d7 ffl = \ cos 26(t)^p- = — ■ The initial value of the phases is chosen as 7j(0) = 0, 
such that 7i(t) = — 72 (^)- For the T-conserving case we find 71,2(0 = 0. The phase 71 (t) 
is determined successively for increasing t from the product of (h(t)\ and In the 

upper panel of Fig. |5]is plotted the phase 71 (t) accumulated by \Ri(t)) when encircling the 
EP twice along the outer loop of the contour in Fig. HI The lower panel shows jx(t) for the 
contour in Fig. HI The start point is the intersection of the loops. The orientation is chosen 
such that the EP is always to the left. The cusps occur where = 0. In each panel the 
triangle marks the start point, the pentagon the point ti, where the EP is encircled once, the 
diamond marks the point t 2 of completion of the second encircling. One sees that in both 
examples 71 (ti) 7^ and that, as predicted in Ref. [lj], the geometric phase 71 (U ) is not real. 



If the EP is encircled twice along the same loop, we obtain 71(^2) — l2Q| within 10 ' 



Thus we measured the 7's up to an accuracy manifestly better than 10 -2 . Accordingly, from 
the value of 71 (£2) given in the caption of Fig. Owe may conclude that, when the loops are 
different 71 (t) does not return to its initial value 20|. When encircling this double loop again 



and again we observe a drift of 7^2 in the complex plane away from the origin. Thus, one 



geometric amplitude increases and the other one decreases. This provides the experimental 
proof for the existence of a reversible geometric pumping process, as predicted in jlj. The 
dependence on the choice of the path signifies that the complex phases 71,2(0 are geometric 
and not topological. In order to obtain the complete phase gathered by \Ri(t)) during the 
encircling of the EP one has to add the topological phase accumulated by \ri(t)). 

In summary, when encircling the EP we obtain geometric factors different from unity 
and the transformation scheme Eq. Q for |ri,2(i))- These results provide an unambiguous 
proof that the EP lies inside both loops of the contour shown in Fig. HI Their precision is 
illustrated by the dense sequence of points in Fig. [51 The accuracy of the measurements at 
and around the EP allows the determination of the size of T violation and of the geometric 
factor along arbitrary contours in the parameter plane. Predictions on geometric amplitudes 
and phases could be confirmed. 
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FIG. 1. Reflection spectra measured at antenna 1 (solid) and 2 (dashed) for s = 1.66 mm, 5 = 
41.50 mm and B = 53 mT. Each antenna couples to predominantly one eigenmode for that choice 
of parameters. We observe a single resonance at antenna 2, a much broader with a shoulder at 
the peak position of the former at antenna 1. Thus one eigenmode is localized in the right part of 
the cavity, whereas the second one penetrates from the left into the right part of the cavity. Inset: 
top view (to scale) of the microwave cavity. In each half of the resonator an antenna, 1 and 2, is 
positioned. A semicircular Teflon disk (gray) is positioned at a distance 5 from the center, s refers 
to the height of the opening (gray bar) between both cavity parts. A ferrite is located at F. 
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FIG. 2. Left panels: Real and imaginary parts of the eigenvalues, fx 2 an d Ti 2, as function of 5 at 
s = Sep = 1-66 mm and B = 53 mT. Around 5 = 5ep = 41.25 mm they are closest. Right panels: 
Modulus and phase of the ratios v\p = |^i,2|e 1,2 of the components of the associated eigenvectors. 
They are also closest at d ~ 41.25 mm. For 5 < 5ep the upper curves correspond to, respectively 
fl, Si, H, $2- 
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FIG. 3. Phase of the ratio vep of the eigenvector components at the EP as a function of the 
magnetization of the ferrite. For the 7~-conserving case the known result [12] of 90° (dashed 
horizontal line) is recovered. The model for the T-violating matrix element iUyi i n t erms of the 
ferromagnetic resonance yields the solid line which agreeably describes the data. The vertical bar 
indicates the range of B where the ferromagnetic resonance is expected. 
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FIG. 4. Differences of the complex eigenvalues fx 2 + iTx,2 ( a ) an d of the ratios Vx,2 = 1^1,2 1^** 1 ' 2 
of the eigenvector components (b) in an area of the parameter plane (s, S) around the EP. The 
darker the colour the smaller is the respective difference. That of fx,2 (blue) and of $12 (orange) 
are small only to the left, those of (orange) and of (green) to the right of the EP. The 
black dots in both panels indicate the chosen contour for the encircling of the EP (see text). 
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FIG. 5. The complex phase 71 (t) for B = 53 mT when encircling the EP twice along the outer loop 
(upper panel), and along the inner and outer loops (lower panel) of the contour shown in Fig. [H 
The triangle marks the start point, the pentagon the point after encircling the EP once (t = t±), 
the diamond that after a second encircling (t = £2)- There, 71(^2) equals (3.123 — i3.474) • 10~ 7 , 
i.e. it is close to zero in the upper case, 0.0931 + f0.0152 in the lower one. 
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